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Abstract 

We present a benchmark in the parameter space of the next-to-minimal supersymmetric standard 
model (NMSSM) that provides for a dramatic multilepton signal and no jets containing 5 or more 
leptons resulting from the cascade decays of the third lightest neutralino, Xs; &iid the lightest 
chargino, Xi^) via light charged sleptons. This is a very clean signal with almost no standard model 
(SM) background. In some cases, a total signal of > 3 leptons -|- jets can be detected at the 
5(T level at the LHC running at ^/s = 7 TeV with approximately 3 fb~^ of data and with less 
than 1 fb~^ when running at ^/s = 14 TeV. In addition, kinematic edges in the invariant mass 
distributions of 2, 3, and 4 leptons are easily detectable with large integrated luminosities (~ 600 
fb~^) which can lead to simple measurements of the mass differences of heavy particles in the decay 
chains, including all combinations of the three lightest neutralinos. 
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I. INTRODUCTION 



Multilepton signals are considered one of the best discovery signals of weak-scale super- 
symmetry at hadron colliders (see, for example, [l-6| and references therein). In the case 
of R-parity conserving supersymmetry, the final state of these decays involve the lightest 
supersymmetric particle (LSP), typically a neutralino, which then results in large missing 
energy associated with the multilepton signal. Many authors have investigated this signal in 
the context of the minimal supersymmetric standard model (MSSM), often working with a 
constrained sets of parameters such as minimal supergravity models 7H9|, gauge-mediated 
supersymmetry-breaking models 10|, or non-universal Higgs masses models A model 
independent approach, in which the particle content of the MSSM is used but all mass 
parameters are taken to be free and independent, was performed in Ref. 12|. 



In these SUSY models, the decay chains are usually initiated by gluino or squark pair 
production, associated production of squarks and gluinos, or the primary production of a 
neutralino and chargino pair. In the latter case, which can be the most relevant when 
squarks and gluinos are heavy, the largest production cross sections are typically from the 
process pp -> — > X2xf when gaugino mass unification is assumed The second 
lightest neutralino state, X25 is favored here over Xi because of its large mixings with the 
Wino and Higgsino gauge eigenstates, which couple to the IV-boson. Both the X2 xf 
then decay down to Xi ^i-nd leptons, neutrinos, and jets via virtual photons, Z and Higgs 
bosons, sleptons, or squarks. This is the nature of the multilepton signal derived from the 
neutralino/chargino sector of the MSSM. 

There exist extensions of the MSSM, however, that contain additional neutralino states 
and therefore allow for the possibility of longer decay chains. The next-to-minimal supersym- 
metric standard model (NMSSM) is one well-known example among various singlet-extended 



models 



ISMlTj. Here, the particle content of the MSSM is extended by one additional gauge- 



singlet, chiral superfield. Its effect on the phenemonology of the model is almost entirely due 
to mixing: the spin-0 singlet S mixes with the Higgs bosons, providing one additional CP- 
even and CP-odd state each, while the spin-1/2 singlino S mixes with the neutral fermionic 
partners of the gauge and Higgs bosons, giving rise to the additional neutralino state. How- 
ever, if the mixing is not large it is possible to have a light, singlinolike Xi while leaving the 
composition of the other neutralinos bearing resemblance to the mixings typically seen in 
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the MSSM. In this case, the dominant sparticle production mode becomes Xsxfi allowing 
for longer decay chains that may not only result in larger trilepton signals but also provide 



for multilepton signals with > 5 leptons. This scenario was first motivated in Ref. 
discussed in the context of the constrained NMSSM (cNMSSM) in Ref. 



18| and 



19|. Similar signals 



20-221. We seek to find a 



resulting from the NMSSM have also been discussed in Refs. 
benchmark in parameter space that results in such a signal and then to demonstrate its 
detectability at the LHC 

In Sec. [Til we discuss the NMSSM, focusing in particular on the neutralino sector of the 
model and how it is different from the usual MSSM case. We then describe how we implement 
this model for the purpose of generating Monte Carlo events in Sec. Illll and discuss the chosen 
benchmark in Sec. IIVI We describe how we model LHC detection and the acceptance cuts 
that we use in Sec. |Vl The backgrounds to our signal are briefly described in Sec. |VT]and the 
detectability of the signal at the LHC and its utility in determining supersymmetric particle 
mass differences is described in Sec. IVIIi A final discussion and summary of our results can 
be found in Sec. IVIIII 

II. NMSSM: AN OVERVIEW 

The NMSSM was introduced as means of alleviating a well-known tension which exists 
in the MSSM, the so-called "/i problem." The MSSM superpotential is given by 

VTmssm = ^-HuHd + uyuQHu - dydQHd - eyeLHd, (1) 

where L and Q are the chiral superfields of the lepton and quark doublets, e, u, and d are the 
chiral superfields of the lepton, up-type quark, and down-type quark singlets, and if„ and Hd 



are the chiral superfields of the two Higgs doublets j23|. Here, /i is the only dimensionful pa- 
rameter. It exists in unbroken supersymmetry and is therefore a supersymmetry-conserving 
parameter that should be of the order of the scale of the complete, unbroken theory. It is 
the only such parameter appearing in the Higgs potential; all other dimensionful parame- 
ters are soft supersymmetry-breaking parameters that should be 0{TeV). For the vacuum 
expectation values (vevs) of the two Higgs states to give sjv^^ + v"^ = vsm = 246 GeV, fj, 
should itself be of the order of the weak scale. This creates a naturalness problem. 

The solution provided by the NMSSM is to generate the fi term dynamically by associ- 
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ating it with the vacuum expectation value of a new field. This is done by removing the n 
term of the MSSM and adding the following two terms to the superpotential: 



+ XSH^H, + -S\ (2) 

Here, 5* is a gauge-singlet, chiral superfield and A and k, are dimensionless parameters of 
order unity ^. Given this superpotential, the effective /i parameter is then given by 

//off = X{S) = A-^. (3) 

All dimensionful parameters of the neutral scalar potential are now 0{TeV) and the singlet 
vev, and therefore /i, becomes naturally of the order of the weak scale. 

In solving the /i problem, the addition of a gauge-singlet, chiral superfield to the MSSM 
alters the mass spectrum of the neutral fields through mixing effects and therefore affects 
collider phenomenology. The spin-0 singlet mixes with and H^, the neutral Higgs boson 
states, to give a total of 3 CP-even and 2 CP-odd Higgs mass eigenstates. The effect on 
the tree level masses, as well as the one-loop corrections derived from the effective Higgs 



potential, are described in detail in Ref. 



25| 



For the multilepton signals, we focus our attention on the neutralino and chargino sector 
of the model. The neutralino mass matrix is given in the (P, W^, H^, H^, S) basis as 



-giVd/2 g,vj2 \ 

g2Vd/'2 -g2Vu/2 

-yUeff -fJ'cSVu/s , (4) 

-/Ucff -fieSVd/s 

-fJ'csVu/s -ficsVd/s V2ks J 

where Mi and M2 are the gaugino mass parameters, and are the up-type and down-type 



f Ml 

M2 

-giVd/2 g2Vd/2 

giVu/2 -g2Vu/2 





Higgs vevs such that ■y^ + i'a 



(246 GeV)^, and tan/3 = Vu/vd- The usual electroweak 



^ The cubic term forbids a continuous Peccei-Quinn symmetry whose spontaneous breaking would introduce 
fine-tuning problems associated with bounds on the non-observation of axions . There remains, however, a 
discrete Z3 symmetry whose spontaneous breaking introduces a cosmological domain-wall problem. This 
is typically circumvented with Planck-suppressed operators that explicitly break the symmetry without 
affecting the weak-scale phenomenology of the theory (see 



241 and references therein) 
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gauge couplings are gi and g2. The upper-left 4x4 is the standard neutralino mass matrix 
of the MSSM, while the outer row and column give the singlino contribution. If the vev s is 
large and V2ks < fics, Mi, M2, the lightest neutralino Xi can be very singlinolike and light, 
with a mass m^o ^ \/2ks. The effects of such a state, when presumed to be the LSP, have 
)een studied in the context of the dark matter relic density [221, 126|, I27| and collider searches 
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2l|. 



The chargino mass matrix is given in the iW"^, ^u/d) basis as 



^ M2 

92Vd 

\ V2 



92Vu \ 

V2 



(5) 



It is left unchanged from the usual MSSM case. With the exception of the modified mixing 
matrices, the couplings of the neutralinos and charginos to the py-boson and sleptons are 
also left unmodified. 

III. MODEL IMPLEMENTATION AND CONSTRAINTS 



To search for a suitable benchmark and to generate our signal, we extended the standard 
MSSM implementation included with MadGraph version 4.4.44 28|. To do this, we included 
the additional neutralino and Higgs boson states and modified all the neutralino and Higgs 
boson couplings to reflect the additional states and mixings and the effects of direct couplings 
to the singlet/singlino states. We independently implement all tree- level sparticle masses 
and mixing matrices. In addition, we include the one-loop effective potential corrections to 
the Higgs boson masses from top / stop loops 251] . All decay widths and branching fractions 
were then calculated using BRIDGE version 2.20 29| and verified with SPheno version 3.0 
jsol . which has hard-coded implementations of the NMSSM branching fractions. 



To generate our signal events, we first calculate the two-body process pp 



X3X1 using our modified MadGraph NMSSM implementation. Using our own Monte Carlo 
code, these events are then decayed down to all possible final states by using the branching 
fractions provided by BRIDGE. These events are then subject to experimental acceptance 
cuts described below in Sec. IVl 

When searching for a benchmark point, we apply a series of constraints by checking each 



parameter set using NMSSMtools version 2.3.2 



3l|. This applies basic collider constraints 



including LEP mass limits, measurements of the anomalous magnetic moment of the muon, 
and upper limits of 6 — 57 as well as theoretical constraints, such as the verification of a 
global minimum of the Higgs potential and the exclusion of points which are found to have 
Landau poles in A, n, ht, or hh when these couplings are run up to Mqut- NMSSMtools also 
includes a relic density constraint of 0.094 < Qh^ < 0.136, calculated using the NMSSM 
implementation of Micromegas 32|. We are interested in points that do not provide too 
much dark matter. Therefore, we only enforce the upper bound Qh'^ < 0.136. In addition to 
the above constraints applied by NMSSMtools, we perform our own independent checks on 
the perturbativity of A, n, and ht, the constraints from muon [g — 2)^, and the LEP limits 
on ZZH couplings. Our implementations are described in Ref. j^^]. 



IV. BENCHMARK 



To search for a benchmark point, we note that we are primarily interested in the produc- 
tion mechanism pp — )■ — )■ xtxf- There are then in principle several decay chains that 
can give rise to multilepton signals with > 3 leptons and no jets. For example, in the case 
of heavy sleptons, the decays may be mediated by real or virtual photons and Z-bosons: 

xl ^ v^*hl ^ i^rv'^*Vi ^ i^i-i'^i'~xl 

where V = A,Z. Another interesting possibility is for a hierarchy of the type M^o > Mj± > 
M^o > Mi± ^ > M^o. The charged sleptons will decay as Z^^^ — )■ l^x^ which then allows for 
the following decay chain: 

xl -> ^ l^l~X2 ^ l^ri'^l'l/L ^ 1^1-1'^1'~X\- 

Similar decays of the chargino are possible if M^± is also larger than the charged slepton 
masses: 

Xt ^ ^il%R ^ ^il^xl ^ ^iI^I'^I'r/l ^ ^il^l'^l''x\- 

Taken together, this decay chain can lead to signals with up to 7 leptons and no jets. In 
order to achieve such a striking signal we seek to generate parameter points that satisfy the 
following criteria: 
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Parameter Scan Ranges 



Lower Limit Upper Limit 



Ml 



— — ivio — — iVl'i 

2 6 
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10 
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2000 GeV 


Ml, Me 


100 


200 GeV 



TABLE I: Parameter scan ranges used to produce the NMSSM benchmark point. The soft scalar 



quark masses are taken to be 2 TeV. The parameter ranges used here are adapted from Ref. 



Model Parameters 


tan/3 




/ieflf K 


A^ 


At Ab Ar Ml Ma Mg Mq Mu Md Ml Me 


7.55 


0.056 488 


199 0.015 


-39.6 


-1170 1886 -143 149 297 891 2000 2000 2000 140 110 



TABLE II: NMSSM model parameters for the benchmark point. The dimensionful parameters 
fieS, A, and M are in GeV. 

• Xs has large Wino and Higgsino components 

• Xi is largely singlino > 0.5) 

• Charged sleptons are light enough that they mediate the neutralino decays 

We proceed with a scan over NMSSM parameter space. We choose our NMSSM-specific 
independent parameter set as: 

s, K,, Af^, Ag. (6) 
We also have the following parameters, which are shared by the MSSM: 
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M(x<() (GeV) 



FIG. 1: The cross section for pp — t- — t- X3X1 ~^ (5 leptons + jets + X) versus the mass of 
the Ughtest neutrahno. The red plus signs represent ah points in the scan, while the blue triangles 
are for those points which have the mass hierarchy M^o,M ± > A/r± > M,i) > A/r± > M,,o and 

A3' Xi 'l a2 'fl Xi 

also |.^]^|^ > 0.5. The yellow circle with the black outline denotes the benchmark point described 
by Table M 



fi,s, tan/3, At, ^, A,, Mi, M^, M3, Mq^, Mu,, M^,, Ml,, Me,, (7) 

where i = 1,2,3 is a generational index. We will suppress this index and assume the 

sfermion mass parameters are the same for each generation. As we are primarily interested 

in neutralinos and light slepton superpartners, we set the squark mass parameters to 2 TeV. 

We also assume gaugino mass unification: Mi = -M2 — -M3. 

2 6 

The first condition from the above list suggests that we take /i > Mi. As described in 
Sec. [ni the second condition often arises when s is very large (greater than several TeV) 
and when \/2ks < min(Mi, M2, /i). Therefore we take k < Mi/ {V2s). Finally, to satisfy 
the last condition we take Ml and Me between 100 and 200 GeV. The remaining parameter 
ranges are adapted from the search for a singlinolike LSP satisfying the relic density used 
in Ref. 26|, where < and Ag > 0. The parameter values and ranges used in our scan 
are defined in Table [B 

We find that large multilepton signals with greater than 5 leptons are fairly generic for 
parameter points in our scan which have the mass hierarchy 

M^o, M^± >Mj±> M^o >Mj±> M^o, (8) 
in addition to a high singlino content of Xi- The scan points are exhibited in Fig. [H which 



Sparticle Mass Spectrum (GeV) Neutralino Composition 



y? • 

Al • 


109 




147 


IJ T D ' 


2020 




B 


w 


Hu Hd S 


A2 • 


129 


■ 


118 




2020 


y? ■ 

Al • 


0.02 < 0.01 0.01 0.01 0.95 




191 


rt ■■ 

i 


114 


h : 


2030 


x^ ■■ 


0.64 


0.03 


0.20 0.09 0.04 


x^ 


206 




150 


h : 


2040 


x^ 


0.33 


0.17 


0.26 0.24 < 0.01 


x^ 


333 


h ■ 


125 


h : 


2010 


x^ 


0.01 


0.01 


0.47 0.51 < 0.01 




173 




125 


h ■■ 


2100 


x^ 


0.01 


0.79 


0.06 0.14 < 0.01 


xt-- 


333 






9 ■■ 


1060 











(a) (b) 



Dominant Leptonic Branching Fractions 
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(c) 



TABLE III: These tables give, for the chosen benchmark, the (a) mass spectrum of the neutrali- 
nos, charginos, sleptons, squarks and ghiino in GeV, (b) neutraUno composition (mixing elements 
squared), and (c) leptonic branching fractions. Here, I = e, fi and Di are the partners of the left- 
chiral neutrino states, u and d represent the scalar partners to the first two generations of up-type 
and down-type quarks, respectively. The states in (c) are ordered according to descending mass. 
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plots a{pp — )■ XsXi" ~^ 5 leptons + jets + X) versus the mass of the lightest 

neutralino. Most points with mass hierarchies described by Eq. [8] and with |.^/^^|^ > 0.5 
have cross sections that are fairly large (~ 1 — 100 fb) and which may be detected at the 
LHC even when including realistic cuts and detector effects. 

We choose a single benchmark point from this scan to be used for a more detailed anal- 
ysis. The values of all model parameters for this benchmark are given in Table HIl The 
resulting mass spectra, neutralino composition, and leptonic branching fractions are given 
in Tables [ml (a)- (c). 



V. DETECTOR SIMULATION AND ACCEPTANCE CUTS 

We organize our signals according to the number of leptons present and enforce a jet 
veto on each event. We do not enforce a Ipx cut or r veto, although these could easily be 
included. Since the signal of interest contains leptons and no jets, we perform our analysis 
on parton-level generated events and use a series of cuts and detector-level effects to roughly 
simulate actual signal detection. Our choices for p^, 77, and Ai? = \J (A^)^ + (A?])^ cuts 
are 



20 GeV for the hardest two leptons (e,/i) 
7 GeV for all other light leptons 
15 GeV for r leptons 
^ 20 GeV for jets 



\ri\ < < 



(9) 



2.4 for electrons 
2.1 for muons 

2.5 for r-leptons and jets 

{0.2 for light leptons 
0.4 for all others 

Detector smearing of the energy of the jets and leptons is modeled as in Eq. [121 



(10) 



(11) 
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y/E/GeY 
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0.03 for jets 
0.007 for leptons 



(12) 



We also include basic acceptance cuts and tagging efficiencies according to Ref. [5| and the 
effect of isolated leptons from heavy quark decay, which occurs with a probability of ~ 1/200 
33|. 

We note that both our signal cross sections and the background cross sections described 
in the next section are calculated at leading order in QCD and do not include the effects 
of showering or initial and final state radiation. As we are ultimately interested in signals 
that do not contain jets at the parton-level and have no colored particles in the primary 
decay chain, the effects of showering would be minimal and most jets from initial state 
radiation would presumably be too soft to be tagged with the above criteria. Initial state 
radiation would, however, have an effect on the resonant VT-boson production cross section 
and kinematics, but as this is a higher-order effect it should not significantly modify our 
results or conclusions. 



VI. BACKGROUND ANALYSIS 

The predominant standard model (SM) backgrounds involve the production and decay 
of weak vector bosons and heavy quarks. These channels are given in Table llVt along with 
their corresponding cross sections at the LHC running at a/s = 7 and 14 TeV when cuts 
and detector effects are taken into account. The parton-level events were calculated using 



ALPGEN version 2.13 3^ while the decays, cuts, and detector effects are later applied. We 
note here that for each background involving an on-shell Z-boson, there is a corresponding 
one with a virtual photon that may also contribute. However, the photon contributions are 
small relative to those from the Z-boson diagrams and we do not include them here. 

The largest backgrounds are naturally in the trilepton channel. These are affected the 
most by the jet veto, as there are nontrivial contributions coming from isolated leptons 
originating from heavy quarks in processes such Zee, Zbb, and tt. The total backgrounds 
for four-lepton signals are also nontrivial but could be greatly reduced with a modest cut on 
^T, as nearly all the signal originates from the leptonic decays of Z-boson pairs. Of note is 
that the backgrounds for final states with greater than 5 leptons are very small, with cross 
sections of (9(10 ab) or less. 

Given these backgrounds with 10 fb~^ of data, estimates of the minimum cross sections 
needed for 3a and 5a signals in various channels, when including acceptance cuts, are given 
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Background Cross Sections (fb) 
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TABLE IV: Multilepton background events for the LHC running at \fs = 7 and 14 TeV. Values 
are given for each lepton multipHcity before and after the jet veto. All values include the cuts and 
detector effects as described in Section |Vl 

in Table |Vl For simplicity, these are calculated in the Gaussian approximation of signal 
significance: 

Significance = ( , = | x -i / / C dt . (13) 

As long as systematic errors on the signal and background cross sections are not large, Eq. [T3] 
should provide a reasonable estimate of the reach. 
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Required Cross Sections (fb) 
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TABLE V: The cross sections for the N lepton + jet signal, after cuts and detector effects, that 
are necessary for 3a evidence and 5a discovery at the LHC running at 7 and 14 TeV with 10 fb~^. 

Signal Cross Sections (fb) 
leptons + jets 
31 Al 51 61 71 
7 TeV 25.6 4.91 2.31 0.09 0.03 
14 TeV 68.7 13.3 6.09 0.29 0.06 

TABLE VI: Multilepton cross sections given by the benchmark NMSSM point for the LHC running 
at both 7 and 14 TeV, after cuts and detector effects are included. 

VII. RESULTS 

A. Signal Rates and Significance 

The cross sections for the signal, at both y/s = 7 and 14 TeV, are given for various lepton 
multiplicities in Table IVII In determining the detectability of the signal at the LHC, we are 
primarily interested in the signals with > 3 leptons + jets and > 5 leptons + jets. The 
total rates for both the signal and the background for these processes at \/s = 7 and 14 TeV 
are given in Table I VI II The approximate luminosities necessary for 3a and 5cr significance 
are also given in this table and are derived from Eq. [131 

At the LHC running at 7 TeV, 3cr evidence can be seen with ^ 1 fb~^ of data while 5a 
discovery is possible with ^ 3 fb~^ in the > 3 leptons + jets channel. Running at 14 TeV, 
discovery at the 5a level is possible with slightly less than 1 fb~^ of data when considering 
the > 3 leptons channel. While the 31 channel is not unique to the NMSSM as it may be a 
discovery channel for the MSSM, the 5/ channel is more exclusive. The LHC may observe 
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TABLE VII: Cross sections for the benchmark signal, and standard model background, in two 
channels (> 3/, > 5/) after accounting for cuts and detector effects for the LHC running at a 
center-of-mass energies of 7 and 14 TeV. Also included are the estimated luminosities required for 
3a and 5a discovery. Here we use I = e, iJ,. 

the 5/ channel with as httle as 10 fb~^ of data running at 7 TeV and 4 fb~^ of data running 
at 14 TeV. While nonobservation of this signal would not allow the NMSSM as a whole to be 
excluded, the mass hierarchy described in Sec. HV] would typically give rise to such a signal. 



B. Kinematic Edge Measurements 

Given such a large multilepton signal, it is then possible to estimate several mass dif- 
ferences from the kinematic edges that will be observable in the invariant mass distribu- 



tions of the leptons when enough integrated luminosity is accumulated |35N42|. Analytic 
expressions for the kinematic edges in the invariant mass distributions of 2, 3, and 4 lep- 
tons resulting from a sequence of two-body decays are given in Appendix |Al For a decay 
chain involving a sequence of two-body decays, each emitting one massive particle and one 
massless lepton, the absolute upper limit of the invariant mass of the leptons is given by 
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I- 1+ I- 1+ 




(c) (d) 

FIG. 2: Primary diagrams contributing to (a) 3 lepton final states (b)/(c) 5 lepton final states and 
(d) 7 lepton final states. 

M^Sixin leptons) = (m/ — mp), where m/ is the mass of the initial mother particle and rriF 
is the mass of the heavy particle resulting from the final two-body decay of the chain. While 
this limit can not always be saturated, it is a fairly good approximation to the true upper 
limit when this mass difference is not too large compared to the masses involved. We can 
therefore interpret these kinematic limits as probes of mass differences of particles at various 
steps in the decays of xt and x^- 

The primary decay chains giving rise to our multilepton signals are illustrated in 
Figs. |2](a)-(d). Each step in the decay of Xs results in the emission of light charged leptons 
via two body decays. Every combination of 2, 3, and 4 leptons in adjacent steps of the decay 
chain therefore offers the potential to measure mass differences by means of kinematic edges. 
The mass differences that may be probed and the true kinematic limits of the associated 
invariant mass distributions are summarized in Table IVIIII It is clear that, in most cases, 
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Process 


Total Mass Difference (GeV) Kinematic Edge (GeV) 


Invariant Mass Distribution 




fi9 A 


oy .o 


2L-0S-SF 


xl^xl 


82.5 


82.5 


4L 


xi ^ li 


72.6 


72.4 


3L 


It ^ 1% 


28.9 


28.0 


2L 


It ^ X? 


38.8 


38.8 


3L 


x^^x? 


20.1 


20.1 


2L-0S-SF 



TABLE VIII: Various subprocesses involved in the decay of as in Fig. [2jd). The total mass 
difference between the initial and final heavy particles of each process are given along with the 
true kinematic upper limit in the associated invariant mass distribution. In this table, I = e, /i. 
The masses and mass differences are given in GeV. The distributions are labeled by the number 
of leptons and, in the case of lepton pairs, whether they are required to be opposite-sign (OS) and 
same- flavor (SF). 

Muiax{n leptons) = (m/ — mp) is a good approximation to the true value. The most useful 
tool is the invariant mass of opposite-sign (OS) and opposite-flavor (OF) lepton pairs, which 
can be used to measure the mass differences M^o — M^o and M^o — M^o. It can also be 
used to measure the mass difference Mr± — Mr±. However, this difference is not limited to 

'l '■r 

opposite-sign, opposite-flavor lepton pairs; the Majorana nature of the neutralinos allows 
their decay into leptons of either charge and flavor and the same is true for the decays of the 
charged sleptons to neutralinos. One of the best ways to measure this edge is then to use 
same-sign dileptons. Beyond the dilepton invariant mass distributions, the invariant mass 
distribution of 3 leptons can be used to measure the differences Mr± — M^o and M^o — Mr± , 
while 4 leptons may be used to measure M^o — M^o . 

The decay of xf offers fewer opportunities to measure differences involving its mass; there 
will always be a neutrino emitted in either the first or second step of the decay chain, resulting 
in a loss of crucial kinematic information. However, when Xi ~^ ^i^ti the subsequent decay 
of the slepton will provide information as described above. 

To demonstrate these effects, we look at various invariant mass distributions of leptons in 
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FIG. 3: Invariant mass distribution of dileptons of (a) opposite-sign, same- flavor leptons and (b) 
leptons of the same sign and for (c) three and (d) four leptons. These distributions are all derived 
from the > 5 leptons + jets channel at the LHC running at t/s =14 TeV with 600 fb~^ of 
integrated luminosity. In all channels, the backgrounds are negligible. The vertical dashed lines 
indicate the kinematic upper limits for various subprocesses which are labeled by the initial and 
final heavy particles involved. In these plots, / = e, /i. 



the > 5 leptons + jet signal with 600 fb~^ of data ^ at the LHC running at ^/s = 14 TeV. 
While offering less signal significance than the case of > 3 leptons due to a smaller overall 
rate, the virtual lack of background makes this signal ideal for kinematic edge searches. We 

^ We choose such a large integrated luminosity to illustrate the kinematic edge effect in the absence of large 
statistical fluctuations in the data. 
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plot the dilepton invariant mass distributions for the > 5 leptons + jet signal in Figs. |3](a)- 
(b) and the three- and four-lepton distributions in Figs.[3](c)-(d). As mentioned above, there 
is essentially no background to these signals and each kinematic edge listed in Table IVIIII 
can be clearly identified. Therefore, not only is the multilepton + jet channel excellent 
for the discovery of this type of NMSSM benchmark, it can easily provide important mass 
information as well. 

VIII. CONCLUSIONS 

We have investigated a scenario in NMSSM parameter space in which a light, singlinolike 
neutralino is the LSP. This allows for extended decay chains from the production pp — )■ 
— >■ xtXi" which can then result in large multilepton signals of high lepton multiplicity. 
In particular, we find that signals with > 5 leptons and jets in the final state can be quite 
large given the mass hierarchy M^o,M^± > M^± > M^o > Mi± > M^o, a singlinolike LSP, 
and heavy squarks. The backgrounds for such processes are virtually negligible at the LHC. 
Therefore, when considering a representative benchmark point and looking at a signal with 
> 3 leptons and jets in the final state, 5a discovery is possible at the LHC running at 
a/s = 7 TeV for approximately 3 fb^^ and it is possible with less than 1 fb~^ of data when 
running at 14 TeV. In addition, the > 5/ + jet channel, which is more unique to this 
benchmark, can be discovered at the 5a level with 10 fb~^ of data running at 7 TeV and 4 
fb~^ of data running at 14 TeV. 

The high multiplicity multilepton signals are also quite useful for measuring mass dif- 
ferences by looking at kinematic mass edges in large amounts of accumulated data. When 
looking at the invariant mass distributions of two, three, and four leptons in the > 5 leptons 
-|- jets signals with 600 fb~^ of data, these kinematic edges are clearly visible, allowing 
for the determination of nearly every mass difference present in the decay of Xs- High mul- 
tiplicity multilepton signals are therefore a very useful tool to use at the LHC should this 
particular scenario of the NMSSM be realized. 

Acknowledgments 

The authors would like to thank the referee for the constructive comments that were 
given and W. Zhu and P. Langacker for their valuable early participation in this study. 
This work was supported in part by the U.S. Department of Energy under Grant Nos. DE- 
IS 



FG02-95ER40896, DE-AC02-06CH11357, DE-FG02-91ER40684 and in part by the Wiscon- 
sin Alumni Research Foundation. 



19 



Appendix A: Kinematic Edge Expressions 



For the decay chain A hB ^ bcC, where each step represents a two-body decay and b 
and c are massless particles, the kinematic upper hmit for M'^ipc) is given by 



{m\ - ml){m% - ml) 



m 



< {niA - mcf 



(Al) 



B 



where mA, and mc are the masses of particles A, 5, and C, respectively 35|. 

For the decay chain A ^ hB ^ bcC — j- bcdD, where b, c, and d are massless particles, 
the kinematic upper limit for M'^ibcd) is given by 
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39|. 



{mA - nin) 

where mA, m^, mc-, and m^ are the masses of particles A, B, C, and D {2], IstI . 

For the decay chain A — t- bB — t- bcC — ?■ bcdD — )■ bcdeE, where b, c, d, and e are massless 
particles, the kinematic upper limit for M'^ibcde) is given by 



Ml^ibcde) 



{m\- ml) {ml -ml) 



m 



B 



{m\ml — mlml){ml — ml) 



2 2 
m^m^j 



{m\ml — mlml){ml — ml) 



mlml 



{m\-ml){ml 



m. 



m 



D 



[mA - rriE) 



iff 
iff 
iff 
iff 

otherwise 



rriA 
rriE 

rriA 
rriE 

ruA 
rriE 

ruA 
rriE 



> 



< 



< 



< 



ml 



ml 



m 



B 



m 



c 



m, 



c 



m 



2 ' 

D 



m 



D 



mi 



where m^, m^, mc, m^, and mE are the masses of particles A, B, C, D, and E 
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